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ABSTRACT 

Since the lower resonant frequencies are of interest in most 
structural problems , a significant reduction of independent variables 
is possible by the use of the normal modes of structural subsystems as 
independent variables. 

This memorandum describes a technique that can be used to 
generate equivalent spring-mass models for the normal modes of a structural 
subsystem when the generalized mass matrix and resonant frequencies are 
available. Where modal truncation is employed, the residual mass matrix 
is used to preserve the correctness of the rigid-body mass properties. 

Applications of the modeling technique and the residual mass 
matrix are discussed. 
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INTRODUCTION 


Since the 3~ver resonant frequencies are of interest in many 
structural problems, a significant reduction of independent variables is 
made possible by the use of the normal modes of structural subsystems as 
independent variables. The representation of the normal modes of the 
subsystems as uncoupled single degree-ot-fre^U >m spring-mass systems simplifies 
the task of combining two structural subsystems which are attached at a 
conmon point and it puts the information in such a form that most structural 
analysis computer programs can be used to evaluate the normal modes of the 
total structure. 

This report uses well-known concepts to develop a technique of obtaining 
an equivalent spring-mass system for each normal mode when the dynamic 
characteristics of the structural subsystem are available as a generalized 
mass matrix and associated resonant frequencies. The dynamic characteristics 
of a continuous subsystem can thus be combined with a discrete system because 
the continuous subsystem can be represented by a set of mutually independent 
single-degree-of -freedom systems. 

A description of the procedure is to renormalize each normal mode such 
that its reactions are represented by those of a corresponding single degree-of- 
freedom equivalent system. Incremental inertia properties must be added in 
the model to represent the rigid-body contribution of ‘•he truncated normal 
modes. To help illustrate the ideas, a simple model is qualitatively described. 

The model is a representation of structure A, a cantilevered beam, attached 
to a structure B as shown in Figure 1. 
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1st Mode 2nd Mode ith Mode 



Fig. 1 INITIAL STRUCTURE Fig. 2 CANTILEVER MODES OF STRUCTURE A 

The equivalent single degree-of-freedom system for each normal mode of 
Structure A (Fig. 2) is normalized and represented such that each reaction, 
R^, is properly simulated when combined with Structure B as shown in Fig. 3. 



Fig. 3 MODELED STRUCTURE 

T’he "rigid mass” in Fig. 3 represents the contribution of the truncated 
modes to the rigid-body mass properties of Structure A. 

The procedure to define equivalent single degree-of-freedom systems is 
developed for the general case wherein six base reactions are represented by 
three orthogonal force components and three orthogonal moment components. 
Appendices 1 and 2 discuss systems with fewer base reactions to aid in the 
explanation of the ideas. 
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ANALYSIS 


The equilibrium Equations 

The equilibrium equations for any discrete or continuous linear, undamped, 
structural subsystem in terms of its rigid-body modes and its arbitrarily- 
normalized characteristic modes are, in matrix form: 


^RR ! ^RKi 


f p ' 

R 


'o 

! o 
! 


N 



1 

■< 

* 

> 4* 

— 

| 

4 



► ^ 

i 


• * 

P 

< r Nj 


0 



l p j 


j 


Where 


( 1 ) 


KJ = mass matrix, rigid-body modes 
T 

KJ = inertial coupling mat i x, rigid-body and normal modes 
[N7IJ = diagonal matrix of generalized masses in normal modes 
IRJ-Ksy = the diagonal stiffness matrix 
| P R | = vector of generalized displacements in rigid-body modes 
fP M l = vector of generalized displacements in normal modes 
(**) = second time derivative of ( ) 

(f} 


generalized force vector 
Eq. (1) may apply to a structural subsystem with inertia relief in some 
rigid-body degrees of freedom, in which case the related N/l^ and 
terms are zero. Here the normal modes are to be regarded as applying to fully 
cantilevered structures. 

Let (L -1,2,?) denote a set of orthogonal reference axes with origin 

at the base point . The vectors, LC- (l - f t 2, V>) , represent translational 


The base point is the point where the model f the structural subsystem under 
discussion is attached to another structural subsystem, effectively at a point. 
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represent rotational 



i*j- 
[<£«] * 


rigid-body transformation matrix 

R. x S matrix of normal-mode disp 1 acements at the point, 
where » and S * number of normal modes . 


The original equilibrium equations can be written in a form also applying 
to a system of single-degree-of-freedom spring-mass elements, ach of which 
represents a normal mode. The following analysis leads to such a system. 
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The Mass Matrix 


For the 



normal mode, 


a renormalization factor is define; 


= (Ml+ H, +- rfjf/ 

(If the numerator is zero, use 

The re normalized elements in the complete matrix are 



It follows that 



( 3 ) 


(4) 


(5) 


If the modes are initially so normalized, then is unity, and subsequent 

application of the criteria [Eqs. (3) and (4)] produces no change. 

The lumped-mass system equivalent to the contim vjs or discrete structural 

subsystem is described in terms of selected normal modes, with the complete 

model represented by the sum of such modal models. 

A mass with magnitude is restrained to move parallel to the vector 

M llf l + K/l i + M, k . The line of action of , is placed such that 

Zhjj "N 

moment arms, , , about the 3^- axes are 

(tt3JN L 






( 6 ) 


* 

Mass terms with a bar represent the original mass terms, and without a bar 
the renormalized mass terms. 
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If = 0 (1 = 1, 2, 3), the location of the rotatory inertia can be 

chosen arbitrarily: however, this inertia ie constrained to rotate about an 

axis parallel to the. vector + Kxj' + ■ 

For the pr nt, a distinction is made between a normal mode of the 

original distributed (D) structure and its equivalent Jumped (L) system by 

D 

use of letter su~ rscripts. The motion of the modal mass, f\/J corresponding 
to a unit rigid-bodv trai slation is 


& = </ M 


P 

NN 


(i = 1, 2, 3) 


For a unit rigid-body rotation about, say, the ^ axis 
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At - !/)%■• - m'/m“ 

^ \ IN // Ntf bib ' Hfi 


Thus, in general. 
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By definition. 
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Moreover, since to 
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( 8 ) 


(9) 


Subsequently, the superscripts, L and D, will be omitted. 
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irom Eqs. (8) and (9) the complete mass matrix for the pi ^ normal-mode 


equivalent may be written as 


1*0 , 
M -j 
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or 







(10b) 
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The Stiffness Matrix 


After renormalization, the stiffness matrix for the 



normal mode is 


where 




0 

0 



I 






( 11 ) 


The original structural subsystem represented by S normal modes can be 

replaced by an equivalent discrete model comprised of S independent single- 

th 

degree-of-freedom spring mass elements. For the Kl mode, the mass is constrained 

to move along a prescribed line of action. Elastic restraint is provided by a 

Z 

spring of stiffness U) H M^. 

The solution can be simplified by eliminating the calculations necessary to 
determine the vector description of the line of action. Moreover, it can be 
made adaptable to structural-analysis computer programs that allow restraints 
in only one coordinate system. 

Accordingly, the displacement system of Eq. (2), which provides for absolute 
base-point motion, is used. 

Let V = base-motion vector 

= modal displacement, in the absolute reference system, of the 

ttl 

point mass representing the n normal mode 
For the single point mass, (Vi,. , 

Nn 

& - 1 

and, for base motion 
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In parallel with Eq. (7), the row matrix 

LM^pJ 


L = - - " > where M 


NfkJ 


is scalar. 



I lolfpj 
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Premultiplying both sides of Eq. (12) by 
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( 12 ) 


gives 


(13a) 


or. 


{p}= [t]{K 


th 


(13b) 


The contribution of the hi normal mode to Eq. (1) is now subjected to a 

th * 

coordinate transformation. For the N normal mode, the new stiffness matrix is 


This came result can be obtained by deriving the force required to give the 
base a unit displacement, successively, in each degree of freedom while holding 
the modal mass, M NN , fixed in the absolute reference system. It can also 
be obtained by requiring null forces for rigid-body motion; i.e., 

0 


r 1 i 



K 1 K 
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LK hb J = -K. w l4vJ 
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(14a) 


Eq. (14a) serves to define the following terms 



(14b) 
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The new mass matrix for the f4 uormal mode is 


r (n) 

M bD 

KR 

M 

! 

titi | 

0 

- 

0 

1 

i 

- 

c 

! 0 
1 



_o 

1 

1 

till - 1 




JPL Techrical Memorandum 33-380 


11 





and the mass matrix 



Eq. (1) becomes transformed to 
where T- 

/<?} - m {f\ 

* See footnote, page 18. 


(lb) 


( 1 */) 

(18) 

(19) 
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The Residual Mass Matrix 


The mass matrix of Eq. (1) is commonly obtained from an initial finite- 
element mathematical model. Its submatrices may be identified as 

( 20 ) 

L J L ' KJ 

i T 


K] - [&] MW 
K„] ■ 

KJ ■ [t] r MW 


( 21 ) 

( 22 ) 


where m is the initial mass matrix of order, and, in general, is non- 

diagoual 

r 1 is the (r X R ) modal matrix for rigid body modes 

L jJ3J 

m is the modal matrix for the normal modes. 

If all f modes are chosen for use in Eq. (1) , M is a square non- 
singular matrix, ard hence can be inverted. 

(*) 

It can be deduced from the form of in Eq. (10) that If all of the 

modes are retained, then 


-I 


[ M J KJ M - [M J 

Indeed, if Eqs. (21) and (22) are substituted in Eq. (23), 

KJ -WlH[i][£] W [t r ] MMM 

r 

* M M [4] (as defined in Eq. 20) 

If, as a matter of engineering judgment, S modes are chosen ( S < r ), 
there is a 1 residual mass matrix," 


(23) 


[^Ws] 


(24) 


Eq. (24) is completely general; its validity does not depend on the renormal- 


ization decribed by Eqs. (3) and (4). 
JPL Technical Memorandum 33-380 
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The residual mass matrix must be added at the base to provide simulation 
of the total rigid-body mass properties. This addition may be done directly 
if the computer program to be used accommodates non-diagonal mass matrices. 
However, two alternate methods are described. 

The first and philosophically most obvious is to choose a sufficient 
number of normal modes to make the residual mass terms negligible. This is 
not necessarily a practical or desirable course. 

An alternative method is to add as many additional spring-mass elements 
as there are rigid-body degrees of freedom; and to assign to each a frequency 


well above the modal frequencies of interest or of validity for the total 
structural system. The [M^] matrix for the added systems is contained in the 


relation 


[Mres] " [^RKl] KJ 


where, for the i, fc normal mode and the L ^ rigid body mode [from Eq.(5)] 

/A- a z 

% = vJr M 4 (2< 

(use A- 4 . . . 6 if £ * 0) 

Eq. (25) may be written as 

KJKJ "KJ 'KJ = KJ (2: 

where the indices, , L , pertain to rigid-body modes and where the row and 
column indices have been included in matrix labels to maintain the 
correspondence to the physical problem. 


The matrix equation (27) can be written as 

KJ - [D]>] 


(28) 


14 


JPL Technical Memorandum 33-380 



Let 


Then 


[D] be the Choleski decomposition of 

K;]- 


The matrix element 



Thus, use of Eqs. (29) and (30) in Eq. (26) gives 


( 29 ) 


(30) 





or 




(31) 


Substitution of Eq. (31) into eier.ents of Eq . (29) gives, with use of Eq. (30), 

'ft* 3 

M;; = / X L V (32) 






(use At- * 4 . . . 6 in Eqs. (31) and (32) if £ * 0) 

Thus the added "high-frequency" modes, of number equal to the rank of 
completely account for the residual mass elements resulting from truncation. 

The generalized masses for these normal modes are obtained from Eq. (31), and the 
"rigid-elastic" coupling matrix is obtained from Eq. (32). 
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APPLICATIONS 

The first documented application of the renormalization concept 

described herein appears in Ref. (1). The only available mathematical 

models of the cantilevered Ranger and Surveyor spacecraft were obtained 

from modal surveys by application of Eq. (14), wherein the O were the 

N 

measured mode shapes associated with physical items of known mass properties. 
In the case of Ranger, only the first cantilever torsion mode was of interest. 
In the case of Surveyor, nine modes having significant coupling with "rigid- 
body roll" were chosen. Since the launch vehicles in both cases were repre- 
sented by discrete torsion-spring and moment-of-inertia elements, it was 
expedient to represent the normal modes of the spacecraft by equivalent 
spring-mass systems. 

The Atlas booster engines and the Centaur main engin s were remodeled 
in a manner to assure only antisymmetric motions. Each engine mode was 
coupled with rigid-tody roll of the axis containing the two gimbal blocks; 
product-of-inertia terms provided the coupling. On renormalization, each 
antisymmetric mode of the engine pair was representable as a single spring- 
mass system, and a residual mass (moment-of-inertia) was attached tc the node 


in the plane of the gimbal blocks to preserve the total rigid-body properties. 
Perhaps the most useful facet of the concepts presented herein is the 


physical feel provided to the analyst by associating the M vector with 

RN 

base reactions and by using the residual mass matrix as a guide in truncation. 


16 


JPL Technical Memorandum 33-380 



As an example, a mathematical model of a spacecraft and its structural 
adapter to the launch vehicle was developed for use with the SA’TS Computer 
Program (Refs. 2 through 5). The mass matrix of this model was of 139th order. 
The first 42 normal modes were computed, along with the and matrices. 
Additionally, the 7 idual mass matrices were com uted for the number of 
consecutive modes from mode 1 through mode 42. Table 1 lists the main- 
diagonal elements, of the residual mass matrices through the 14th 

mode, and lists, also, the A associated with each added mode, in the 
first mode, the generalized mass associated with the (translational) 
coordinate (i.e., M 2 2 ) accounts for 82% of the rigid-body mass; and the 
generalized mass associated with the x^ (rotational) coordinate (i.e., M^) 
accounts for 72% of the rigid-body moment-of-inertia about the x^ axis 
(see Fig. 4) . 

The first ten modes contain little effective mass in the longitudinal 
(x^) direction; the 11th and 12th modes account for 85% of the rigid-body 
mass effective in the x^ direction. 

TheZ^K^ of Table l permit qualitative description of the character 
of a me without any detailed knowledge of mode shape. Those modes for 
which all A are small are termed "local modes," in which some component 
or appendage is the principal contributor. 

This type of information can be of use to the analyst who is undertaking 
a dynamic loads analysis of the entire space vehicle. If the obj active is 
to obtain a preliminary assessment of loads on the launch vehicle, economies 
can be effected by deleting "local modes," with placement of the appropriate 
residual masses at the base of the spacecraft. 
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During a modal survey of the Mariner Mars : 71 Development Test Model, 
an on-site computer terminal was used to make an orthogonality check of 
measured modes afte T completion of the survey of the second and each subse- 
quent mode. The residual mass matrix was also computed. Inspection of this 
matrix served as a guide in shaker placement for excitation of a new mode. 

It also served as a practical indicator of the number of modes to be surveyed. 

CONCLUSION 

The concepts presented herein are intended to aid the physical under- 
standing of the dynamic influence of a complex substr w rc at its point of 
attachment to another substructure. There may be instai in which it is 
convenient, tor one reason or another, to adopt the equivalent spring-mass 

approach through the described renormalization process. 

* 

The concept of the residual mass matrix, which is not dependent on 
any particular modal normalization, is, nonetheless, a byproduct of the 
modeling concept. It has proved to be very usefu 1 in the exercise of 
engineering judgments relating to modal truncation in analysis and to the 
requisite completeness of a modal vibration survey. 


*Footnote: 

After the completion of the draft of this memorandum, it came to the 
authors' attention that Schwendler Slid M&cNeal, in Ref. 11, define and use 
"residual flexibility matrices." Moreover, in writings not in the open 
literature, MacNeal has used the term "residual mass matrix" with exactly 
the same def ition as gi\ an rein. 
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Table 1. Residual Mass Matrix Elements Applying to Consecutive Numbera of Normal Modes of a Spacecraft Model 
Having an Initial Dynamic Matrix of 139tn Order. 
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APPENDIX 1 


111 ostrative Case of a Single Base-Reaction Component 


Consider the longitudinal modes of a "fixed-free" uniform bar of length 
and weight per unit lengthy r . From Ref. (6), it can be shown (using 
Timoshenko's notation) that the circular frequency of the ith normal mode is 



where 





( 1 . 1 ) 


The modal amplitude at a distance % from the root is 

X = D, (i *{,%£■■■) a-2) 

where is, in general, an arbitrary normalization factor. 

Here it will be shown that, for each mode, a value of C^‘ exists to give 
base reaction equivalence between the continuous system and a simple spring-mass 
model . 

The generalized mass of the ith mode is 


M„ = f 
U 0 

where 1 Q &!• (1*3) leads to 

M. • 


(1.3) 


( 1 . 4 ) 


The axial load in the ith mode at station ?C is 


p i f i / 


(1.5) 
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where X r is an arbitrary rigid-body displacement. 

The integral in Eq. (1 C) is recognized as the "rigid-elastic" coupling term 
in the partitioned matrix 


wherein 




[rvl r , : 

i 

M r f 

X: 
^ 1 
J 

Rl ;i . 


- vL 

] yL f Md% 


With 



1 . 0 , 


m -yn 

i i rr 

P L - f! X 


(1.7) 


Use of Eq. (1.2) ir the Integral of Eq. (1.6) gives 



IjF' 


( 1 . 8 ) 


To replace the ith mode of the continuous system with a simple system as 
sketched, it is necessary only to specify a modal renormalization that equates 
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the inertia force to the force acting on the base: 



i. e. 



Thus, equating (1.4) and (1.8) gives 



and tne mass matrix becomes 



j 



"m„!_ 

K : 



• - • M 

1 Vr7 




M tr ! 
. 1 


Kt 

$ 

I 



' M„ n 


It is convenient to rewrite Eq. (1-11) as 



%YY1 


(c -1,2,3,...) 


(1.9) 


( 1 . 10 ) 


( 1 . 11 ) 


( 1 . 12 ) 
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For the infinite set of elastic modes 


L - cO 


r K; - 


L - OO 

Snl y [ 

-Tf 


From Ref. (7) 


/ 



(1.13) 



(1.14) 


From Ref. (8), the sum of the infinite series 


!~h 






(1.15) 


Adding t.ht«e series gives 



(1.16) 


(1.17) 


In most practical engineering situations, use of a relatively few of the 
lower modes suffices. Thus, in order to keep the rigid-body mass properties 
correct, a residual mass is added to the base. 
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The equivalent mathematical model of the " 
in Figure 1.1, which also shows the decrease in 
number of modes selected. 
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APPENDIX i 


Illustrative Case of Two Base-Reaction Components 

Consider the bending modes of a cantilevered uniform beam, for which 
there will be transverse shear and bending moment to be reacted by the base. 

References 9 and 10 deal with characteristic functions of uniform 
beams with various boundary conditions. The modal normalization has been 


chosen such that, for the nth mode. 


/ k (x ) d* 


( 2 . 1 ) 


where is the beam length. Letu^ be the mass per unit length. Then the 
generalized mass in the nth mode is 


M nn = j 4> n (*}dx 


( 2 . 2 ) 


The shear at the beam root is 



if - 


A i 


(2.3) 


where 



is an arbitrary rigid-body translation. 
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In this context , 

the integral of Eq. (2.3) 

is the "translation-elastic 

" coupling 

term, M-j-^ in 

v ' ,ass matrix 





' Kr 

m tr m ti 

Re - 

• M 

T4/ 
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^RR 

Rz • • 

• R* 


m - 

M 
,Y, i T 

M ie 



(2.4) 


M ir 

9 

M 2 |? 
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R* 




. Kr 
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The bending moment at 

the root is 
a 





0? -- 
r 

*1 (zfi„ 

zjaz 




O 

* 0 


(2.5) 


where ~ ^ the transver&c displacement of an element of the beam at 

station^ for an arbitrary rotation, , about the root. Thus, the integral 

of Eq. (2.5) is the "rotation-elastic” coupling term 

J 

Rn * | < 2 - 6 > 

0 
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and 


Now 


Since the normalization of th° 


rigid-body 


moues is optional, let 



A 


be unity. Then 






> 


(2.7) 


( 2 . 8 ) 


(2.9) 


From x. tegral 1 of Ref. 10 



From Integral 25 of Ref. 10 





( 2 . 10 ) 


( 2 . 11 ) 
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Now, to represent the nth normal mode by a simple spring-mass system. 


first choose a renormalization factor 



Then the new matrix elements are 



(The shear at the base is now "matched. ") 



(2.13) 


The base moment may be matched by placing the lumped mass at station 



M 




Mhn 



(2.14) 
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In the model, this placement may be done in at least two ways: 



M 


nn 


K no = CO* M 


n 'nd 





3 

n 


It can be shown that 


Ilf?! 




(2.15) 


In r 





(2.16) 


N-»°° 




I 



(' . 17) 
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For a finite number of modes, M , a residual mass must be attached to the 


base. 


M 


Rets. 


IV'tt ~ 


H 

y \a 

^ ' ' tin 


( 2 . 18 ) 


By placing this mass at 



(2.19) 


the static-moment equivalence is preserved. 

In general, the rigid-body moment-of-inertia equivalence may require the 
addition of a centroidal moment-of-inertia. 


( T) = VI - 

v O'kK 


, n 'Res, Res. 


or 


(Ij 


o 'R€$ 


M 




£ + M«(^ 


(2 




m 


The [M] and [ Mj matrices for the first five normal modes of a ’.aaiform 
cantilever beam are presented below: 


. 20 ) 
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[ m ] - m 


M- 


yyi 


1 i 

2 ! 

f I 

.78299 

.56883/ 

.43394 .25443 .18190 

.090767/ .032416/ .016542/ 

.14147 

.OlOOO?/ 

~T“ 

1 

1 

1 

1 

1 

1 

1 

0 0 0 

0 

1 

1 

1 

1 

1 


10 0 

0 

1 

1 

1 


1 0 

0 

(SYM) 



i 

i 

i 


1 

0 

i 

i 

i 

i 

i 



1 

/ I 

.61307 

.18830 .064735 .033088 

.020014 

2, j 




/ ! 
3 j 

.44539/ 

.039387/ . 0082476/. 0030089/ .0014157/ 

i 

i 

i 

i 

i 

i 

.61307 

0 0 0 

0 

i 

i 

i 

i 

i 


.18830 0 0 

0 

« 

i 

i 

i 


.064735 0 

0 

i 

i 

1 

(SYM) 



1 

i 

i 

i 


.033088 

0 

i 

i 

i 

i 

i 



.020014 

l 





( 2.2 


( 2.2 
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Table 2.1 gives a description of residual mass requirements as a function 
of the number, N, of uniform-beam elastic modes chosen. It illustrates the 
approach of the residual mass to zero as N increases. 

A model representation is given in Fig. 2.1. 


N 

— 

M 

RES 

X 

RES 

(I _ ) 

RES 

M 

TT 

i 

M 

RR 

1 

.3869 

.1411 

.00616 

2 

.1986 

.0767 

.00103 

3 

.1339 

.0521 

.00031 

4 

.1008 

.0394 

.00010 

5 

.0808 

.0316 

.00004 


TABLE 2.1 RESIDUAL luASS DESCRIPTIONS AS A FUNCTION OF 
THE NUMBER OF ELASTIC MODES CHOSEN. 
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The stiffness matrix associated with this system is 



From numerical values of 



listed in Ref. 


9 


the stiffness matrix 


for the first five modes is 


[K] = K 




.6131 


7.395 


19.33 


39.12 


64.66 


where 



Z 

aJ l m 
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Application of the coordinate transformation of Eqs . (14^ and (15) leads to 



i 

i 

i 


and 


/%] = m 



J 


I 


A 






131.7 12.66/ ; 

-.6131 

-7.395 

-19.93 

-39.12 

-64.67 : 

1.618/ j 

i 

i 

-.4454 / 

-1.547/" 

-2.540./ 

-3.556/ 

4 

-4.574/; 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

j 

.6131 

7.395 

19.94 


1 

i 

i 

i 

j 

(SYM) 



39.12 


i 

j 

! 

1 





64.67 


.0808 .00255/ 

.00255/ .OOOlo/ 

0 

0 

.6131 

.1883 | 

.06474 | 


.03309 


.02001 
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As a farther illustration, assume that only the first three normal modes of 
the uniform beam are to be chosen, and that the truncated modes are to be represented 
by two additional spring mass systems of arbitrarily high frequency. 

From Eq. (28), the Chcleski decomposition matrix is written as the upper 
triangular matrix 



(T: translation) 
(R: rotation) 


The residual mass matrix, determined from values listed in either Eq. (2.21) 
or (2.22) by use of Eq„ (24), is 




a 

IZ 


.13389 

.00697 

C2-2_, 

1 

3 


.00697 

.00047 


Then, by Eq. (28) 
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D 4f 



CL 

II 

^"12. 


-m 


O 

1 

04 

df 



^22- 


JPL Technical Memorandum 33-380 


37 



frcn? which 





D D --&-•/ 

U 4 f-R/ Hi 



> 






^2f 

#"11 


From Eq. (32), and the above. 



From Eq. (31), 


2 






2 . 


Dks e 

^22 


Z 

du 

£L[i 
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In matrix form, for the two added degrees of freedom. 



Using these M;; , M:; , and matrices to replace the last 

oo qL 

two rows and columns of Eq. (2.22), and making the coordinate transformation 


indicated by Eqs. (15) and (14) leads to 
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and, with the original frequency ratios of modes 4 and 5 arbitrarily retained for 
the two residual-mass modes , 



186.3 

12.77./ 

-.6131 

-7.395 

-19.93 -158.31 0 



1.729// 

. 

-.4454 / 

. 

-I.547/ 

-2.54o/ -8.24l/ -.3463 / 

II 

V 

i 1 

(SYM) 


.6131 

7.395 

19.93 

158.31 






“ . 

.3463 •£ 


L 
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